Shift of Dirac points and strain induced pseudo-magnetic field in graphene 
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We propose that the strain induced effective pseudo-magnetic field in graphene can also be ex- 
plained by a curl movement of the Dirac points, if the Dirac points can be regarded as a slowly 
varying function of position. We also prove that the Dirac points must be confined within two 
triangles, each one has 1/8 the area of the Brillouin z( 
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The discovery of graphene, a monolayer carbon atom 
sheet [H , and the development of experimental technique 
to manipulate this two-dimensional(2D) material have ig- 
nited intense interest in this system [2|-l5|. One of the 
most attractive characters of graphene is that its low 
energy excitation satisfies a massless 2D Dirac equation 
@, and the chemical potential crosses its Dirac points (or 
Fermi points) in neutral graphene. These special char- 
acters lead to many unusual properties and new phe- 
nomena 0, •> such as the anomalous integer quantum 
Hall effect (QHE) [8, 9]. Recently, experiments have con- 
firmed another remarkable effect that mechanical strain 
can induce a very strong effective pseudo-magnetic field, 
leading to a pseudo-QHE, which can be observed in zero 
magnetic field [Io|, Gj] • In this paper we propose that the 
strain induced effective vector potential can also be ex- 
plained by shift 5K(x) of the Dirac points K(x), its effec- 
tive pseudo-magnetic field is in proportion to V x K(x), 
only if the Dirac points K(x) can be regarded as a slowly 
varying function of position, and the Fermi velocity is 
generalized to a tensor [12]. We also prove that the Dirac 
points can not be arbitrarily moved, they must be con- 
fined within two triangles, each one has 1/8 the area of 
the Brillouin zone(BZ). 

Firstly, consider a tight-binding Hamiltonian describ- 
ing a uniformly deformed honeycomb lattice with three 
different nearest-neighbor hopping energies ^1,^2,^3 na- 
na = 



<iajfc> 



h.c, 



(1) 



where Cj b (c\ a ) are annihilation(creation) operators, i(j) 
are position vectors of unit cells, a(b) denote two inequiv- 
alent atoms in a unit cell, ti a ,j6 is the electronic hopping 
energy from the jth unit cell b atom to ith unit cell a 
atom. Suppose that the deformed lattice remains invari- 
ant under spatial translation, i.e., U a jb only depends on 
i — j, but the three nearest-neighbor hopping energies 
£1,2,3 ma y be different owing to anisotropy of strains, as 
shown in FigHJ The hopping parameters can be written 
as some 2x2 matrixes t(i — j), whose elements are de- 
fined by [t(i — j)] a ,b = Uajb- For this nearest-neighbor 




FIG. 1: Unit cell and hopping parameters for deformed 
graphene. 



tight-binding Hamiltonian, the non-vanishing hopping 
matrixes are 



t(0) 



£1 
ti 



,t(ai) 



t 2 




,t(a 2 ) 



t 3 




(2) 



and t(— ai) = t^(ai), t(— a 2 ) = tt(a 2 ). By Fourier trans- 
formation 



~7^^2 c Ka(b) exp(ik-j) 



Mb) 



with N a normalization constant, the Hamiltonian (pQ) 
can be cast into the form 



h aa {k) ftoft(k) 
h ba (k) h bb (k) 



Ck,a 



where h aa (k) = h bb (k) = 0, h ab (k) = ft£ a (k), and 
h ba (k) =ti+t 2 exp(zk • ai) + £3 exp(zk • a 2 ), 



(3) 



(4) 



with ai,a2 the lattice unit vectors. The energy bands 
obtained by diagonalizing this Hamiltonian are|l6j 



£7±(k) = ±|ti+t 2 (k) + t3(k)|, 



(5) 



where t 2 (k) = t 2 e ikai , ? 3 (k) = t 3 e ika2 , the plus sign 
corresponds to the upper (tt) and minus to the lower (jr*) 
band respectively. From Eq.((5]) we notice that if K is a 
zero point of h ba (K)^ i.e., 



ti+f 2 (K)+f 3 (K) = 0, 



(6) 



2 



then £+(k) and £L(k) will meet at K, i.e., E+(K) = 
£L(K) = 0, this K is known as the Dirac point. The 
Hamiltonian ([3]) can be expanded up to a linear order in 
p = k — K in a neighborhood of point K 



h aa (k) h ab (k) 
h ba (k) h bb (k) 



a* -p 

dp o 



(7) 



where \i, v — 1,2 denote two components of a 2D vector 
and a sum over the repeated indices \i, v is implied, a 
is a complex vector with Re(df) = (^11,^12), Im(a) = 
(^21, ^22), cr 1 ' 2 are Pauli matrixes acting on the sublattice 
degree of freedom, tensor v^ v represents the anisotropy 
of the dispersion near the Dirac points, it only occurs 
noticeable departure from vpS^ in a strongly deformed 
graphene[12|. However, after this modification the strain 
induced effective vector potential will acquire a direct 
physical meaning. For a graphene under nonuniform but 
slowly varying strain, *^(x) and hence the Dirac point 
K(x) as well as v^^x) can be regarded as some smooth 
functions of position x, the local linearized Hamiltonian 
r^(k u -K u (x.)) on the RHS of Eq.Q can be cast into 



tV(xK(p"-«T(x)), 



where p - £K(x) = k - K(x), 5K(x) = K(x) - K f 
with Kj the corresponding Dirac point in strain- free 
graphene. Unlike the usual explanation of the strain in- 
duced gauge field in graphene 0, QjJ, where the effective 
vector potential is an auxiliary quantity and describes 
the mixed effects of both anisotropy of v^ v and the shift 
of Dirac point, here the vector potential only represents 
the relative translation of the Dirac points, (e/c)A(x) = 
5K(x), its pseudo- magnetic field B(x) = (c/e)V x K(x), 
and the physical effects are mainly determined by the 
pseudo-magnetic flux through a loop (c/e) § L K(x) • dx. 
In the following sections we shall discuss the properties 
of K(x), and illustrate how a curl field K(x) is induced 
by a strain. 

From Eq.© we know that the vectors representing *i, 
*2(K), ?3(K) in the complex plane can form a directed 
triangle for a Dirac point K, as illustrated in Figj2ji. 
According to the triangle inequality, we have the follow- 
ing necessary and sufficient conditions for the existence 
of the Dirac points [13j: 



*i + * 2 > * 3 , * 2 + * 3 > *i, * 3 + *i > * 2 . 



(9) 



These conditions define a pyramidal domain in the 
(^1,^2^3) space, shown in Fig|2]3. If *i,* 2 ,*3 satisfy in- 
equalities (j9j), then there exists two directed triangles 
with the same edges *i,*2,*3 but different possible ori- 
entations, which determine two angles 61,62 satisfying 
*i + t 2 e %ei + ^e* 6 * 2 = 0, where 6\, 6 2 are given by the law 
of cosine 

f 2 _|_ +1 j-2 



6±i — ±[tt — arccos (— ^ 2 



6± 2 = ± [ arccos ( — 



2*1*2 
*3 ~ to \ 



(10) 



2*1*3 



7T . 




FIG. 2: (color online), (a) Zero points of /i& a (k) determine 
two directed triangles with edges 1 1,2,3 in the complex plane. 
For a given *i and a fixed direction of *2, the arguments of 
*3 must satisfy conditions ([T5]) to ensure *2,3 > 0. (b) Dirac 
points exist if U satisfy inequalities (|9]) , which describe a pyra- 
midal domain in (*i, *2, *3) space, if (*i, *2, *3) goes beyond this 
domain, an energy gap will be opened. 



(8) Thus the Dirac points K can be determined by letting 



exp(iK • ai) = exp(z#i), exp(zK • 8l 2 ) = exp(z^2), (11) 



so we have 



K 



2tt 



(#ibi+# 2 b 2 )+Ko, 



(12) 



with bi,b2 the reciprocal lattice vectors defined by 
• bj = 27r5ij, and Ko = nbi + mb2 with n,m are 
arbitrary integers. Notice that if *i + i 2 + *3 = 0, then 



*i 
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*3 =0, this implies that there exists two Dirac 
points K(x) and — K(x). However, if (*i,*2,*3) exactly 
locates on the boundary surface of the pyramid, e.g., 
*i = *2 +*3, then the two triangles will mutually coincide 
and * 2 = *5, *3 = *3(see FigEK), hence K(x) and -K(x) 
become equivalent, and a = i(t 2 a.i-\-ts8i 2 ) becomes a pure 
imaginary vector, so the Fermi velocity in the directions 
perpendicular to a vanishes (Fig J3b and l3H) [13. HEI. fljjl . If 
(*i, *2, *3) goes beyond the domain defined by Eq.([9]), e.g., 
*i > *2 + *3, Eq.((6]) will have no any root, an energy gap 
with magnitude E g = 2(*i — 1 2 — ts) will occur at the cor- 
responding points K± = d=l/2(bi +b2)(see Figl3b)[l9|. 
and the effective Hamiltonian (j8]) must be further modi- 
fied by adding a mass term and some second order terms. 

Another important property is the range of K(x). We 
shall prove that the Dirac points must be confined within 
some special regions of the BZ. To this end, notice that 
if a Dirac point K = (l/27r)(0i, 6 2 ) is given, then its 
associated *i,2,3 can also be determined up to an arbitrary 
factor, except six special cases of 61,62 = 0, ±7r(see Fig. 
E^l). If 61,62 ^ 0, ±7r. According to the law of sines we 
have 



t_2 

*i 



sin #2 



sin(<9i -6 2 y 



*i 



smt/i 



sin t/ 2 



(13) 
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FIG. 3: (color online) (a) Energy band when two Dirac points 
are very close, where t\ = 2.8, £2,3 — 1-45, (b), (d) £1 = 
2.8, £2,3 = 1.4, two Dirac points are equivalent (superposed), 
(c) £1 = 2.8, £2,3 = 1-35, an energy gap occurs. 

or 

(£i,£ 2 ,£ 3 ) oc (sin(0 2 - 0i), - sin0 2 , sin0i). (14) 

For the six special cases we have: if (61,62) = ±(71", 7r), 
h = £ 2 ±fe; if (0i, 2 ) = ±(tt,0), £2 = ti+£ 3 ; if (0i,0 2 ) = 
±(0, 7r), £3 = £1 + £2- From Eq.JTH) we can find that the 
0i , 02 must satisfy some constrain conditions to guarantee 
£1,2,3 > 0, as illustrated in Fig|2ji. For an arbitrary £2 and 
a fixed 6\ (direction of £2), £3 must point in a direction 
between the directions of —£2 and negative real axis, i.e., 
argument 0i , 02 must satisfy 

01+7T<0 2 <7T, 01 e (-7T,0), 
-7T<0 2 <01-7T, 0iG(O,7r). ^ ' 

These two inequalities respectively determine the range 
of K(x) and — K(x). They describe two open trian- 
gles AMM[M" and AM'MiMg in reciprocal space, as 
shown in FigHJ each one has 1/8 the area of a unit cell of 
the reciprocal space(the parallelogram M" M'^M^M^), 
and each Dirac point is confined within a triangle, so, 
the Dirac points K and — K can meet (become equivalent) 
only at the vertexes of AMM[M" and AM'MiMg. The 
remaining hexagon(blue in Fig0]) is a forbidden region 
for the Dirac points. This confinement also limits the 
order of magnitude of V x K(x), i.e., the strain induced 
pseudo-magnetic field. In order to show the underlying 
regularity, here we have ignored the variations of bi,b2 
with the deformation of lattice, and simply sketch all 
K = (fci,fc2) in the same affine frame. After translating 
to the first BZ of graphene, AMM[M" and AM'M\M'{ 
are equivalent to a ringlike region consists of six triangles 
AMKM[, AM1K2M2 1 AM" K\M' , etc. 




FIG. 4: (color inline) Rang of the Dirac points consists of six 
triangles in first the BZ, or AMM[M" and AM'MiMg . 

In order to illustrate how a non- vanishing V x K is in- 
duced by strain, we only need to analyze three ideal cases, 
in which only one U is slightly changed, U — >• £0 + 5U, 
while the other two tj^ remain constant, tj = £& = £0, 
which can also be roughly regarded as that the bond 
is elongated(or compressed) while the other two bonds 
Cj,Cfc and their directions remain fixed(see FigJSb). No- 
tice that the Dirac points only depend on the relative 
proportions of £i,£2,^3, so, as an equivalent case, we 
can always assume that £1 remains constant and only 
£2, £3 are variables. Moreover, in these equivalent cases 
the £2 and £3 can be determined by the end of the vec- 
tor £1 ± £2 = —£3, denoted by P in FigJSti. So, we 
can represent the variation of the Dirac points by the 
shift of the point P. To this end, we have to determine 
the corresponding P of the three classes of characteris- 
tic points in the range of the Dirac points: (1) if (or 
if') etc. (see FigH]), according to Eq. (fT3]h Dirac points 
locate at these two points only if £1 = £2 = £3, their 
corresponding P is located at if (or if') in FigJSJi; (2) 
critical points M(M 1 ) = (±1/2,0), M'(M[) = (0,±l/2), 
and M"(M!{) = (±1/2, ±1/2), in these cases there ex- 
ist only one Dirac point since the points K and — K are 
equivalent, their corresponding (£i,£2,fe) are located on 
the boundary of the pyramidal domain, while their cor- 
responding P are located at the real axis in Figj5^; (3) 
O, O', O" etc., their corresponding P are the centers of 
two circles and the infinite limit points of the straight 
line KK' in FigJS^i, which respectively correspond to 
the limits of £2 — » 0, £3 —> and £1 — » (equivalent to 
t2 = ts — » 00). Now we analyze the shifts of the Dirac 
points in the three ideal situations. (1) £1, £2 remain con- 
stant and £2 = £1, only £3 is variable, the trajectory of the 
corresponding P is a circle with radius £1 and centered 
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1 — >a: 




FIG. 5: (a) The trajectories of P in three ideal cases, (b) A 
schematic diagram of the shift of Dirac point S~K and its curl, 
5~K perpendicular to c*, if bond a is slightly elongated. 



at the point (ti,0) in Figj5t, so the arguments (61,62) 
satisfy 



9i-2(9 2 -2tt = 0, 6>iG(-tt,0), 
h -2<9 2 + 2tt = 0, 0i G (0,tt). 



They describe line segments M'0\ and M{0 in FiglU (2) 
fe(= £1) remain constant while £2 is variable, the trajec- 
tory of corresponding P is another circle with radius t\ 
centered at the origin, its associated (61,62) satisfy 



9 2 -2(9i + 2tt = 0, 62 e (-tt,0), 
9 2 -2(9i-2tt = 0, 6> 2 e(0,7r), 



(17) 



which describe MO' and M\0[ in Fig|4j (3) t\ remains 
constant while £2, £3 are variable but £2 = £3 (or vice versa, 
ti is variable, £3) remain constant), the trajectory 

of P is straight line KK' , (61,62) satisfy 



9 2 =0, -<|0 1 |<7T, 



(18) 



which describe M"0" and M3 O" in FigH Summariz- 
ing Eqs. (fTBj) (fTT|) (fT8j) and comparing with FigHJ we ob- 
serve that if a band, e.g., ci is slightly elongated (or 
compressed) along its direction, ci —> (1 + S)ci, while 
the other two bonds 02,03 remain fixed, then t± will 
be slightly changed while ^(= £3) remain constant, the 
Dirac point K will be slightly moved in the direction 
perpendicular to ci, i.e., SK y 7^ 0(see FigHJ K moves 
towards O" if t\ decreases, towards M" if t\ increases), 
K'(= — K) is moved in the opposite direction. Thus, if 
the elongation of ci is slowly varying in the x-direction, 
i.e, dti/dx 7^ 0, then dK y /dx ^ 0, the other two cases 
are similar. So, a nonuniform strain as schematically 
shown in FigJSb can induce a curl field K(x), VxK/0. 
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